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Abstract 

The O(n) spin model in two dimensions may equivalently be formu- 
lated as a loop model, and then mapped to a height model which is con- 
jectured to flow under the renormalization group to a conformal field 
theory (CFT). At the critical point, the order n terms in the partition 
function and correlation functions describe single self-avoiding loops. We 
investigate the ensemble of these self-avoiding loops using twist operators, 
which count loops which wind non-trivially around them with a factor 
— 1. These turn out to have level two null states and hence their corre- 
lators satisfy a set of partial differential equations. We show that partly- 
connected parts of the four point function count the expected number of 
loops which separate one pair of points from the other pair, and find an 
explicit expression for this. We argue that the differential equation sat- 
isfied by these expectation values should have an interpretation in terms 
of a stochastic(Schramm)-Loewner evolution (SLE K ) process with re = 6. 
The two point function in a simply connected domain satisfies a closely 
related set of equations. We solve these and hence calculate the expected 
number of single loops which separate both points from the boundary. 



1 Introduction 

In a recent paper, Werner [Q has shown there exists a measure on simple loops 
on any Riemann surface which has the property of conformal restriction. This 
is to say that if D' C D are any two subdomains of the manifold, the measures 
on loops in D' obtained by (a) restriction of the measure on loops in D to those 
in D', and (b) conformally mapping D —> D' ', are the same. Moreover this 
measure is unique up to multiplication by a constant. We refer to these loops 
throughout the paper as self-avoiding loops and to the mass of any subset under 
Werner's measure as the /i-mass. 

One may also consider the set of self-avoiding polygons on some regular lattice 
embedded in the manifold. The total number of such polygons of length I is 
known to grow as /J where fj, is lattice dependent. The measure which weights 
each polygon with a factor x l c (where x c = n^ 1 ) has the restriction property, 
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and is commonly conjectured also to be conformally invariant in the limit of 
vanishing lattice spacing, which means that it should give a particular example of 
Werner's measure. In this paper we assume this to be true, and hence conjecture 
values of the /^-mass of certain loop subsets, using (non-rigorous) Coulomb gas 
and CFT methods applied to a generalised model, the O(n) model. 

The O(n) model encompasses a large group of physically interesting models, 
including the Ising spin model, the percolation problem and self-avoiding loops. 
The theory may be written as a spin model with nearest neighbour interactions, 
or alternatively in the loop gas picture, in which the states of the model are 
non-intersecting loops constructed on the edges of the spin model. Each loop 
is weighted with a factor nx l , where x is related to the reduced coupling in the 
model and I is the length of the loop. In the Ising model, these loops are the 
cluster boundaries of a spin model defined on the dual lattice. At the critical 
point of the model, x = x c (n), it is conjectured to flow under the renormalization 
group to a Gaussian free field 0131 , which is a well known example of a conformal 
field theory (CFT). One of us 0] showed the existence of operators in this theory 
whose two point correlation functions count loops around one of two points with 
a weight n' rather than n, hence giving information about the distribution of 
loops. The choice n' = —n is of particular interest; these operators are then 
called twist operators and are at (r,s) = (1,2) in the Kac classification 0. 
Hence, they have null states at level two. The two point correlation function is 
then equivalent to the expectation value of (— l) w in the loop ensemble, where 
N is the number of intersections of loops with a defect line, a continuous curve 
connecting the two points. The twist operators may therefore be thought of as 
a source and a sink for this defect line. 

Self-avoiding loops are described by the loop gas picture of the O(n) model with 
n — > 0. The partition function and correlation functions are determined (to 
order n 1 ) by graphs with just a single loop. The two point correlation function 
of twist operators therefore leads to an analytic expression for the number of 
single loops which separate the locations of the operators, weighted by x c as 
in figure n This number is logarithmically divergent in the continuum limit, 
however, due to the contribution from vanishingly small loops around each of 
the points. 

In this paper, we consider the four point correlation function of these twist 
operators on the Riemann sphere. In the loop gas picture, we may think of 
each operator as again being a source (or sink) of a defect line. Hence there is 
a defect line running between each pair of points. We argue that the choice of 
paths for the defect lines is unimportant and that they may run between any 
two distinct pairs of points. The four point correlation function can then be 
shown to be the expectation value of (— 1) where TV is now the total number 
of crossings of loops with defect lines. Particular semi-connected parts of this 
four point function yield the expected numbers of weighted loops (the /i-mass 
of loops) which wind around two of the four locations of the operators, as in 
figure |2J CFT may be used to derive the form of the four point function, 
since the null states of the operators imply that the correlation functions obey 
a set of Bclavin, Polyakov and Zamalodchikov (BPZ) type partial differential 
equations [Hj- We solve these equations and hence find analytic expressions for 
the ^t-mass of loops which separate one pair of points from the other. For loops 
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Figure 1: Two examples of single loops winding around one of two operators 
for the theory on the Riemann sphere. The stars mark the locations of the 
operators and the dashed line indicates that the curve wraps around the back 
of the sphere. Two point function is the sum of such loops, weighted by x c to 
the power of their length. Notice that the sum contains loops with vanishingly 
small perimeter in the continuum limit. 

which separate (zi,Z2) from (2:3,2:4), for example, we obtain the expression 



where r\ = (z\ — 22X^3 — z±)(z\ — 23) 1 (2:2 — 24) 1 is the cross-ratio of the points. 

These expected numbers of weighted loops are finite in the continuum limit 
because there is no contribution from vanishingly small loops, and are invariant 
under conformal transformations. We show that the non-leading behaviour of 
these expressions as rj — > reveals the derivative of the central charge with 
respect to n at n = 0. This is consistent with the interpretation |Jj of the stress 
tensor as the spin-two component of the relative probability that a loop passes 
between two points in the limit Z12 — > 0. 

We use similar arguments for the O(n) model in a simply connected domain. 
Such a domain may always be mapped via a conformal transformation to the 
upper half plane with the real axis as the boundary. The two point function 
in this domain satisfies the same partial differential equations as the four point 
function in the bulk, with the locations of the two operators and the reflections 
of these points in the real axis being the locations of the four operators in the 
bulk theory. We show that the connected two point function of twist operators in 
the n — ■> limit counts the mass of loops around the locations of both operators, 
and find the explicit expression 
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where r\ is now the cross ratio of the two points together with their conformal 
images in the boundary. 




Figure 2: Two examples of single loops winding around the locations of two 
of four operators for the theory on the Ricmann sphere. A particular semi- 
connected four point correlation function counts the number of such loops, 
weighted by x c to the power of their length. Note that, in contrast to the 
two point correlation function above, there are no vanishingly small loops con- 
tributing to the sum. 

The layout of this paper is as follows. In the next section, we recall the ar- 
guments leading to the Coulomb gas picture of the O(n) model. In section |3 
we discuss twist operators in the O(n) model and derive the form of the two 
point correlation function of these operators. The small n limit of the theory is 
explored and the correlation function in this limit is shown to depend (to order 
n 1 ) only on the configurations of a single loop. In section 0] we discuss the four 
point correlation function of twist operators. Because the twist operators have 
null states at level two, this correlation function satisfies a set of partial dif- 
ferential equations, which we solve analytically. The interpretation of the four 
point function in terms of a pair of defect lines is explained, followed by the 
example of the four point function of twist operators as the four point function 
of spin operators in the Ising model. Section l4~Tl is devoted to the small n limit 
of the four point function and contains a derivation of the mass of loops winding 
around two of the four points. We then explain how these numbers, as func- 
tions of the positions of the four points, hold a key to measuring the effective 
central charge of the O(n) model as n — > 0. In section [B] we apply the theory 
of twist operators to the O(n) model in a simply connected domain and show 
that the mass of loops around two points in such a domain is invariant under 
conformal transformations and finite in the limit of vanishing lattice spacing. 
Finally, in section [7] we show that the mass of loops around two of four points 
on the Riemann sphere can be thought of in terms of an SLEg process. 

2 The Coulomb Gas 

Let us start with the following partition function for the O(n) model 

Z (n) =TrU(l+xs(r i )-s(r j )). (1) 

m 



4 



The s(rj) are n-component spins on a lattice {r^} and the product is over pairs of 
nearest neighbours. The product in the partition function can be expanded into 
a sum of 2 N terms, where N is the number of nearest neighbours. Each term 
can be represented by a graph of open and closed loops in the following way: 
the neighbouring sites, i and j, are joined together by a line if the xs(ri) ■ s(rj) 
term was chosen in the expansion of equation or left disconnected if the 
number 1 was chosen instead. The trace over spins may then be done for each 
graph individually. Since Trs(ri) odd = because of the symmetry under the 
transformation s(rj) — > — s(r,), all graphs containing an odd number of s(r^) 
make no contribution to the partition function. These are the graphs with open 
loops, hence the partition function will only contain contributions from graphs 
with closed loops. Let us consider the honeycomb lattice for simplicity, then 
there may be 0,1,2 or 3 powers of a given s(rt). The trace over spins of s(ri)° 
will contribute a factor of 1, whereas, from Tr s a (r)sb(r) = 5 a b, instances of 
s(r 2 ) 2 lead to a factor of n for each closed loop in the graph. The partition 
function is therefore equivalent to 

Z 0(n) = (2) 
G 

where I is the total length of all the loops in a given graph, v is the number of 
closed loops and the trace is now over all graphs of closed, non-intersecting loops. 
This form of the partition function may be used for general values of n, whereas 
equation Q is applicable only to positive, integer n. There exists a critical point 
in the theory, at x = x c , where the mean loop length diverges and the model is 
supposed to become conformally invariant. In order to formulate a field theory, 
the non-local factors of n must be made local. This can be done by assigning 
orientations to the loops (either clockwise or anticlockwise) and inserting local 
factors of e lx at each vertex where the curve turns to the right and e~ %x where 
the curve turns to the left. Summing over the two possible orientations for 
each loop leads to a contribution of e &%x + e~ 6lx from each closed loop on the 
honeycomb lattice, where closed loops have a difference in the number of left 
and right turns of six with the sign depending on the orientation of the loop. 
In order to obtain the desired factor of n for each closed loop, x must therefore 
be chosen such that 

e 6*x + e -6* x ^2cos6x = n. (3) 

This may now be transformed into a height model. The height variables are 
taken to be integer multiples of 7r and are assigned to sites on the dual lattice, 
such that the loops are contours of the landscape. Crossing a loop running from 
left to right leads to a decrease in the height of it, whilst crossing a loop running 
from right to left leads to an increase in the height of tt. A given configuration of 
the heights corresponds to a unique graph of orientated loops. The assumption 
of the Coulomb gas is that this height model flows under the renormalisation 
group (RG) into a free field theory with action 

S[h(r)} = £ J (dh(r)fd 2 r, 

for some g{n). There is an additional subtlety regarding topology, which may be 
seen by considering the model defined on the cylinder. On a cylinder of circum- 
ference L, the loops wrapping around the circumference are counted incorrectly 
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(for ri 5^ 2). Loops wrapping around have the same number of left turns as 
right turns. These loops are therefore not counted correctly. The situation is 
remedied by placing a vertex operator (also known as an electric charge) e %%xh l^ 
at one end of the cylinder and an operator e~ l6xh ^ at the other end. Let these 
points be ±w/2. If there is a single loop wrapping around the cylinder, then 
there will be a height difference of tt between the ends, the sign of which depends 
on the orientation of the loop. Summing over the two possibilities for the ori- 
entation then leads to the required contribution of e l6x + e~ l6x . These charges 
at the ends of the cylinder lead to a modification of the partition function to 

^Coulomb Gas = Zff c ( e ~ l ~' l ( _t0 / 2) e l — h ( w / 2 ) ) ffc 



Zffc 
Zffc 
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where the abbreviation ffc signifies that the above expectation values and par- 
tition functions are for the free field theory on the cylinder and where we have 
used w 3> L. The free energy per unit length on the cylinder for this Coulomb 
gas (CG) partition function may be calculated as 

Fcg _ hajZcG) _ TTc 
if w 6L 

This is the form of the free energy on the cylinder obtained via the Schwartzian 
derivative from a theory on the plane with central charge c jS] given by 
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where the first term comes from the known behaviour of Zff c pj. All that 
remains is to fix the value of g, which may be obtained from the following 
argument: adding a term —A J cos(2 h)d 2 r to the action should not affect the 
critical behaviour, since the height variables were defined to be integer multiples 
of 7r and cos(2m7r) = 1 for all integer m. Hence, this term should be marginal 
under renormalisation group flow and therefore must have scaling dimension 2. 
This determines |B] 

ff = l-6xA- ( 5 ) 



3 Twist operators 

In the loop (Coulomb) gas picture of the O(n) model, loops wrapping around the 
cylinder may be counted with a weight n! different from n by placing additional 
charges e ±eih ^ x ~ x " 7r at the ends of the cylinder, with \' chosen according 
to equation Q- The scaling dimension of these operators may be obtained 
from their two point correlation function in the loop gas ensemble. Placing the 
charges at ±s/2, 

( e -^( 1 (-»/2) e -.^i(- s /2) e ,^ft( s /2) e ,&k(„/2)) tac Md ^ 

-™ (e -^M-/2) e ^M-/2) )frcc _ ficld cyL ■ 
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These free field expectation values may be calculated explicitly since they are 
Gaussian moments. Then, in the limit s>i, the two point function is seen to 

4tt 2 

This expression is of the general form of a two point function of operators on 
the cylinder with scaling dimension [H] 

x(n,n')-6 2 (x' 2 -X 2 )/27r 2 .g. (6) 

A particular choice of interest is n' = —n, for which the scaling dimension is 
x(n,—n) = 3/2g — 1; this is the scaling dimension of (f>i_2 operators with level 
two null states in the Kac classification: 

= (rg - s) 2 (g l) 2 

These correspond, in string theory language, to operators which insert orbifold 
points corresponding to the global symmetry s — *■ — s of the hamiltonian and 
are known as twist operators. They will be the focus of this paper. Their 
correlation functions may be considered in geometries other than the cylinder. 
For example, on the Riemann sphere, the two point function of such operators 
is fixed by scale invariance to be 



(</>(zi,Zl)</>(z2,Z2))cFT = — — , (7) 

a 

where we have inserted explicit factors of the lattice spacing, a, so as to make 
4> dimcnsionlcss. The two point correlation function in the loop gas may also 
be calculated on the sphere using the ensemble of equation J5J). Each loop 
separating the points (z\,~z\) and (za>^) will be counted with weight — n rather 
than n, hence graphs in G will be weighted with an additional factor of (— 1) 
to the power of the number of loops separating the two points. Graphs in G 
with an odd number of loops separating point (zi,z\) from point (22,22) will 
be weighted by an additional factor of (— l) odd = — 1 } whilst those with an even 
number of loops will be unaffected, since (— l) ovcn = 1. These two different 
types of graphs can be separated in the sum over graphs as follows 



2-3/fl(n) 



-'loop gas 



(4>(zi,z 1 )(f>{z 2 ,z 2 ))i oop sas = \ ^2 h l ) x l nV + 

G od d(zi,Z2) G ovon (zi,z 2 ) 

(8) 

where G dd(<Zii z 2) is the set of graphs of closed, non-intersecting loops with an 
odd number of loops separating point (zi,zT) and (22,^2) and the G e ven(^ij ^2) 
is the set with an even number of loops separating the two points. This two 
point correlation function has a natural interpretation in terms of a defect line 
joining points (z\, ~z\) and ( 22,^2)- If is the number of times loops cross the 
defect line in a given graph, then equation (jSJ) may be rewritten as 

{4>{Zl,Zlj4>{Z2,~Z5))\ oop gas = ^ {~l) Nl2 x\n m J Zi 00 p gas 

G 

= ((-l) Wl2 ) 1 oopgas. (9) 
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Figure 3: A graph on the Riemann sphere belonging to {G dd} with one loop 
(solid line) and two choices of path for the defect line (dashed line). In (a) 
N 12 = 1 so that (-l) Nl2 = -1. In (b) N 12 = 3 so that (-l)^ 12 = -1 also. Any 
choice of path for the defect line will lead to (— l)^ 12 = — 1 since N12 is always 
odd. 



The defect line can take any path between the two points because, for a given 
loop configuration, (-T)^ 12 will be the same for all paths. This can be best 
understood with reference to figure^ Combining equations Q and ©, we find 



((-l^hoop; 



2-3/ S (n) 



(10) 



3.1 The small n limit 

Self- avoiding loops correspond to g = 3/2, which is the dilute phase of the 
n — > limit of the 0(n) model. The results of section [3] are summarized by 
equation (|1C)|> . which for n«l may be expanded in powers of n. The order n 1 
term in the loop gas expansion (the left hand side of the equation) comes from 
graphs with a single loop. By equating this with the order n 1 term of the right 
hand side, we shall demonstrate a property of the sum over graphs with a single 
loop. 

The left hand side of equation (|10|l for n«l becomes 
(£(-l)^4n° + ) ( 4n° + £ xW) 1 + 0(n 2 ) , (11) 

Go Gi Go Gi 

where Go is the set of all graphs with no loops and G\ is the set of all graphs 
with one loop, hence their respective powers of n° and n . Go contains a single 
graph with no edges, hence the first term in both brackets is equal to 1. The set 
G\ contains graphs with only a single loop. This set is composed of two subsets 
with no overlap, the set Gi jZl | Z2 with a single loop separating the point (zi,z\) 
from (z2,Z2) and the set G\ tZlZ2 \ which contains no loops separating the points. 
It should be noted that a single loop surrounding the point (2:1,21) is classified 
in the same group as a loop around (^2,^2) since on the Riemann sphere one 
can be continuously deformed into the other. Both have an odd number of 
intersections with a defect line between the two points, hence (— l)^ 12 = — 1. 
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For all graphs belonging to G l ziZ2 |, it may be seen that (— l)^ 12 = 1. This is 
shown in figure 0] 




Figure 4: An example of a graph belonging to Gi jZlZ2 | and one belonging to 
Ci, 2l | 2 , for a choice of (z\,Zi) and (22,^2) on the Riemann sphere. Notice that 
a loop in Gi )Zl i Z2 can be thought of as being either around (zi,zi) or around 
(z2,~z~2)- Similarly, a loop in Gi jZ1Z2 i can be thought of as being around both 
points or neither. 

Equation (jl 1|> can therefore be rewritten as 

(1+ e (-i) x y+ e E E + 

G l,*l> 2 G l,-l*2l G l.-ll*2 G L-122l 

= l-2n E < + 0{n 2 ). (12) 

Equations (JBJ, (JSJ and J5J may be combined and then expanded in powers of 
n to find the following form for the scaling dimension of the twist operators in 
the small n limit: 



x(n) = -1 + 3/2g(n) = n/3n + 0{n 2 ) . 
Therefore, the expansion of the right hand side of equation (|1U|I is 



z\ - z 2 



-2n/37r 



+ 0(n 2 ) = 1 



2n 
3tt 



In 



z\ - z 2 



Q(n 2 



(13) 



Hence, from comparing the coefficients of n 1 in equations (|12(l and (|13|1 . we may 
conclude that 



E 

G l,z, U 



3tt 



Zl - 22 



(14) 



1 1=2 



In words, for a given pair of points (z\,zi) and (z 2 , ~z 2 ), the number of weighted 
loops which separate the two points, weighted by x l c , diverges as (l/37r) In \l/a\ 
in the limit of the lattice spacing tending to zero. The cause of this is the 
diverging contribution from vanishingly small loops in the continuum limit. The 
factor x l c may be thought of as the measure, hence this weighted number of loops 
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is equivalent to the expected number of loops around one point. It may be noted 
that the coefficient l/3ir differs from that seen on the annulus with shrinking 
internal radius. In the case of the annulus with vanishingly small modulus |10| 
a coefficient of 1 /6-7T is seen, because there are only diverging contributions from 
small loops around one point, as opposed to the plane where small loops around 
both points contribute. 



4 The four point correlation function of twist 
operators from conformal field theory 

Twist operators were introduced in section [21 as the operators responsible for 
counting loops with weight — n rather than n in the loop gas picture of the 
0(n) model. Their scaling dimension was calculated and it was seen that it 
corresponds to the scaling dimension of operators with null states at level two. 
Conformal field theory may be used to derive a set of partial differential equa- 
tions satisfied by the correlation functions of such null state operators. In the 
case of the four point function, these differential equations may be solved ana- 
lytically, as will be seen in this section. To make contact with the language of 
Schramm (stochastic)-Loewner Evolution (SLE), g is dropped in favour of the 
parameter k used in SLE The two are related by the formula k = A/g. 

It was shown in section [3] that the twist operators have a null state at level 
two U3- This null state is itself a highest weight state; it is annihilated by all 
raising operators, L n with n > 0. This is not the same as saying that the 
null state decouples from all other states in the theory, as it docs in unitary 
theories. That said, it can be shown by a modular invariance argument that 
physical results are seen only if the null state decouples on the torus. This is 
because, barring unforeseen cancellations, the nondecoupling of the null states 
would lead to a density of states at high energies corresponding to c = 1 rather 
than c < 1 22- We shall therefore set the null state to zero in this theory and 
use the resulting partial differential equations for the correlation functions of 
the twist operators. In the complex plane, the four point function satisfies the 
following set of partial differential equations (for j = 1, 2, 3,4) 



d% - J^2~ ( — J 2 ^ 2 f 1 — (</>(zi> z l)</ , ( 2: 2, z 2 )<j>{z 3 , Z 3 )0(z 4 , Zi )) = 



(15) 



hi h d 

<3f-- (—^—\2 ~ —1*— ^{ z i^)^ z 2^)<P{ z z,zi)<P{ z i^)) = 0. 



4 (z,[ Zj) %% 

J 



(16) 



In terms of Zij = Zi — Zj, the cross ratios are defined as 

^12^34 _ Z12 Z 3i 

v = n = = ■ 

Z13Z24 Z13 Z24 

These cross ratios are invariant under global conformal transformations. The 
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partial differential equations above have the solutions 

Z12Z34Z23ZU 

(17) 

where G(rj, 77, re) is 

^,r7 )K ) = | 2J F 1 (l-|,2-^;2-|;r ? )| 2 

+B(K)Kl-,)| a *-| a F 1 (5 > ^-l;5 ; ^ J 

/12 = 3re/16 — 1/2 and h% = k/2 — lis the value of the scaling dimension of an 
operator with a null state at level three. The form of B{k) can be derived by 
requiring consistency under the permutation of the labels (zi,~Zi) and is found 
to be PIUI) 



B(k) 



[r(i - f ) 2 r(f ) 2 - r(2 - f ) 2 r(§ - i) 2 ]r(^ - i) 2 



r(f) 2 r(f - i) 2 r(i 



k\2 



The form of A(k) is dependent on the normalisation of the two point function 
since, as rj — > 0, the four point function becomes 

(cl)(z 1 ,z^)(j)(z2,z^)(f>(z 3 ,z^)(j)(z4,z^))cFT -* A(k)\ 12 2 34 | 4,12 • 



4.1 Interpretation of the four point correlation function 

Recall from section [21 that the correlation function of two <j> operators has an 
interpretation in terms of a defect line in the loop gas ensemble. The operators 
are the source and sink of the defect line, which can take any path between the 
two points. Each graph in the sum is weighted by an additional factor (— l)^ 12 
where N12 is the number of intersections of loops in that graph with the defect 
line. There is a similar interpretation of the four point function in the loop gas 
ensemble, but now there are two defect lines. The four points are split into two 
pairs and defect lines run between the points in each pair. For the choice of 
pairing (21,27) with (22,22) and (23,23) with (2:4,24), the four point function 
can be seen to correspond to the following expectation value in the loop gas 
picture 

{<f>(zi, Zl)<f>(z 2 , Z2~)(f>(Z3, Zs)4>{Zi, Zl))loop gas = ((-l)" 12 ( - 1 ) ) loop gas , (18) 

where N\2 is the number of crossings of loops across a defect line from (zi,zT) 
to (22, 22) and N34 is similarly defined. There are three ways in which the pairs 
may be chosen, but each choice leads to the same set of weights for the graphs 
in the partition function. The reason for this may be seen in figure El Thus, 
the four point function may also be considered to be the expectation value 
in the loop gas ensemble of (-l) iVl3 (-l) iV ' 24 or (-l) Nl4 (-l) N23 . Comparing 
equations (|17J) and (1181) leads to the main result of this section 

((-l)^(-l)^) loopgas = Ma 2 ih > A(Knri ^ Kh (19) 

212^34^23214 

where rj, k) is defined below equation (1171) in section^ 
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Figure 5: (a) shows the choice of defect lines described in the text. As the defect 
lines may take any path between their endpoints, they may be re-routed to pass 
infinitcsimally close to each other, as in (b). Since a loop crossing a defect line 
leads to an additional weight (—1) in the loop gas partition function, a pair of 
defect lines leads to a factor (— l) 2 = 1 and hence makes no contribution to the 
weight of any graph. Hence, the defect lines shown in (b) are equivalent to the 
choice in (c). A similar construction may be used to show the equivalence to a 
pair of defect lines drawn from (21,21) to (22,22) and from (23,23) to (24,24). 
Note that this feature is a result of the choice n' = —n in section 



4.2 Example: the Ising model 



Boundaries between Ising spin clusters correspond to the loops of the 0(n) 
model at n = 1, or K = 3. The twist operators in this case are equivalent to 
magnetisation operators because the parity of the number of loops separating 
two spins depends on whether they are parallel or anti-parallel. That is to say 
that 



((-1) 



N 12 



I loop gas 



Hi , 



21 - Z 2 
(s(2i,2T)s(2 2 , 22"))l s ing 



Similarly, the four point correlation function of the twist operators is the four 
point function of the magnetisation operators. The following function is the 
result of substituting k = 3 into the results of section 01 



A(3) 

(0(2i, zr)0(2 2 ,22")0(2 3 , 23)0(24, 24)) = 



2132240 



2l22 3 42232l4 



1/4 



(|l + ^| + |l-0y|). 

(20) 



Equation (|20|l is the well known four point correlation of spin operators in the 
Ising model at criticality 



5 The small n limit of the four point function 

In section we examined the small n limit of the two point function of twist 
operators, equation I|1U|) . We found an exact expression for the /Lt-mass of loops 
surrounding one of the two points. This result is summarised by equation (|14|l . 
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A similar expansion to order n 1 of the four point correlation function (equa- 
tion I19f) yields information about the configurations of a single self-avoiding 
loop around four points. 

From the equations in section the following small n expansions may be de- 
duced 

A n + 0(n2) 

9n 



±n + 0(n 2 ) 

— n + 0(n 2 ) . 



K = 


3 + 




n 


h 2 = 






Gtt 




1 


h 3 = 


3 + 


B = 





The expansion of A(n) around n = is A(n) = 1 + gn + 0{n 2 ) where g is a 
constant independent of n, since the correlation function should tend to 1 as 
n — ► 0. The Taylor series expansion of the right hand side of equation (|19|l is 
then 

(</>(zi, ^T)^(Z2, ^2)0(^3, z£)(f>(Z4, 2I))CFT = 

- In 1772:232:140 ~ 2 | - ^^773^2(1,1, -; 2,-; rj) + ^3^2(1,1, g5 2 ^>V) 



3V3r(i) 2 r(|) 2 " v 1 v 3' '3 



0(n 2 ). (21) 



The left hand side of equation i|19|) may also be expanded in powers of n. The 
graphs up to order n 1 in G may be split into the graph with no loops (n°) and 
eight distinct sets of graphs with a single loop, as will be shown in the next 
section. 



5.1 The configurations of a single loop around four points 

In this section, we expand the left hand side of equation (|19fl for n <C 1 to order 



n 



The coefficient of n 1 may then be equated with the result for expansion of the 
right hand side, described in the previous section and leading to equation (j21(l . 
Recall that the partition function for the loop gas expansion is equation (0 . The 
sum may be decomposed into a term of order n° coming from the graph with 
no loops (hence the zeroth power of n) and graphs with a single loop weighted 
by nr. Other graphs will contribute terms of order n 2 or smaller. The graphs 
with a single loop on the Riemann sphere may be further split up into eight 
distinct subsets, which are shown in figure El Each configuration has a unique 
set of (-1)^ where i ^ j and i, j G {1,2,3,4}. 

Let L be defined as 

/ 1 \7V 12 { 1 \ W34 T1 #loops™.total length of loops 

T = I ( ( _-\\N 3 4.\ _ 2^G I 1 7 I 1 7 n x c 

U > 1 > /loop gas- ^ n #loop Sa ,total length of loops 
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Figure 6: The eight distinct configurations of a loop around a given set of four 
points on the Riemann sphere. 
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The numerator of L can be calculated as 



Numerator = l + n[ E (1)(1)4+ E + E C 1 )^ 1 )^ 

f~1 s~i f~1 

z l Z2 Z 3 Z 4 " Z 1 Z 2 Z 3\ Z A Kj ' z 1^2 z a\ z 3 

+ E E (-!)(-!) 

J'-l 

y - r z 1 z 3 z 4: \z2 ~1 3 2 1*3*4 "*1*3 1*2*4 

+ E E (-l)(l)^c] +0(n 2 ). 

The denominator takes the same form as the expression above, but without the 
factors of (1) and (—1). For ease of notation, define 

W ZlZ2Z3 z 4 \ ^ ] x c 1 

and similar expressions for the other graphs in figurc[B] The elements of { Wi} are 
sums over all loops of a single configuration type i, weighted by x c to the power 
of their length. They therefore represent the /i-masses of loops of configuration 
i. It may then be seen that 



L = 1 - 2n 



W Z1Z2Z3]Z4 + W ZlZ2Zi[z3 + W ZlZ3Zi[z2 + W Zllz2Z3Zi + 0(n 2 ) (22) 



The two point functions Cr^} = (4>(zi, Zi)(f>(zj, Zj)) = ((—l) NiJ ) and products 
of pairs of two point functions can be expanded in terms of the {Wi} also: 

C 12 C 3 4 = 1 - 2n(W z + W Z + W Z 

IZ2Z3 \ Z4 

+ W Z1]Z2Z3Z4 + 2W ZlZi]z2Z3 + 2W ZlZ3]z2Zi ) + 0(n 2 ) (23) 
C 13 C 24 = 1 - 2n(W z 

1Z3Z4\Z2 

+ W Z 

IZ2 24 I Z3 

+ w z 

+ w z 

1 \Z2Z3Z4 

+ 2W Z 

\Z4\Z2Zs 

+ 2W Z 

1Z2 \Z3 Z4 

)+0(n 2 ) (24) 
C14C23 = 1 - 2n(W ZlZ3Z4]z2 + W Z1Z2Z4]Z3 + W Z1Z2Z3 \ Z4 

+ W Zllz2Z3Z4 + 2W Z1Z2{Z3Z4 + 2W Z1Z3{Z2Z4 ) + 0(n 2 ) . (25) 

Thus, it can be seen from equations (|2*2*|) . f2l^l . f2*4^l and (J2EJ that a subset {W c } 
of the { Wi} can be expressed in terms of partly-connected four point functions: 

W Z1Z4 \ Z2Z3 =-^(L~ C 12 C 3i - C 13 C 2i + C U C 23 ) + 0(n) (26) 
W Z1Z3 \ Z2Z4 =^(L~ C 12 C 3i + C 13 C 2i - C 1A C 23 ) + 0(n) (27) 



W Z1Z2 \ Z3Z4 = ^( L + Ci 2 C 3i - C 13 C 2i - C U C 23 ) + 0{n) . 



These {W^} correspond to graphs with a loop winding around two of the four 
points and hence do not include contributions from vanishingly small loops in 
the continuum limit. By comparison with the small n expansions of L (section|3J) 
and C{jj} (section 13.11 with the same choice for the normalisation of the twist 
operator as in the four point function) from conformal field theory, these weights 
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are the following functions of 77, 77 

-1 



Gtt 

where 



W ZlZi{z2Z3 = —ki\l-ri\+q{r),rj)+0{n) (29) 
W ZlZ3 \ Z2Zi =q{n,r))+0{n) (30) 
W ZlZ2 \ Z3Zl =^ z \n\n\+q{r),r}) + 0{n), (31) 



q{V,v) = 24^^3^2(1, 1, -; 2, -; r?) + rj 3 F 2 (l, 1, 3! 2, -; rj) 

Ml-v)\ 2/3 \2FAl,ir,)f 



3V3r(i) 2 r(|) 2 " v 1 v 3' '3' 

Note that the W c above arc finite in the continuum limit of vanishing lattice 
spacing. In fact, the expressions in equations Q29|l . I |3()}l and (|31|l arc independent 
of a. They are finite and non-zero in the limit n — ► and are invariant under 
conformal transformations, being functions only of the cross ratios. These ex- 
pressions for the /j-masses {W c } constitute one of the main results of this paper. 



5.2 The central charge 

It is a standard result of conformal field theory that, given the explicit form 
of any four-point function, the central charge c may be determined. This is 
because the operator product expansion (OPE) of any scalar primary operator 
cj> with itself takes the form 

4>(z u zi)<l>{z2, z 2 ) = \z 12 \- 4h (l + • • • + ™T( Zl ) + • • • J (32) 

where T is the holomorphic component of the stress tensor. The coefficient 
follows from consideration of the limit z±2 — *• in the three point function 
(T(z)4>(zi)(f>(z 2 )) oc 2h and the two-point function {T(z)T(z 1 )) = (c/2)(z - 

If this result is applied to the explicit expression l|17fl . wc find that the coefficient 
of rj 2 is 



h 2 {2h 2 + 1) 



2fa(l-f)(2-2f) , (l-f)(2-f)(2-f)(3-f) 



2-f 2(2 - f )(2 - f ) 

and so obtain the known expression for the central charge in terms of k 

_ (3k-8)(6-k) 
C ~ 2k ' 



(33) 



This result is of interest in the limit n — > for the light it sheds on the inter- 
pretation of the stress tensor T in terms of an observable of a random curve 
given by Doyon, Riva and Cardy [Jj. In that paper it was shown that, for any 
measure on simple random curves which satisfies conformal restriction, one may 
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identify T(z) as being proportional to the spin-2 angular Fourier component of 
the probability that the curve intersects a small line segment of length e centred 
at z. While in that paper the focus was on curves which connect two points 
on the boundary of a simply connected domain, described in the case when 
conformal restriction holds by SLE 8 / 3 , the theorem should equally well apply 
to Werner's measure on self-avoiding loops, if suitably re-interpreted in terms 
of //-masses rather than probabilities. However our results in the present paper 
do not directly apply to the intersection with a line segment, but rather to the 
event that the loop passes (or does not pass) between pairs of points (which, 
however, may be taken to mark the ends of the line segment). Nevertheless, one 
might expect these to differ only by a constant of proportionality, and indeed 
our results in this paper support this, and suggest that the two-point function 
of the object T introduced in Ref. [7] is indeed given by the central charge for 
small n as expected. 

Indeed, let us consider the quantity W ZlZ4 \ Z2Z3 as given in equation 1291 This is 
the /i-mass of loops which separate (2:1,2:4) from (22,23). In the limit 212 — > 0, 
writing Z\<x = ee 1 ® 12 , we can define 

V( Zl ;z 3 ,z 4 ) = ^ lime- 2 J d6 12 e- 2 ^W ZlZ4lz2Z3 , (34) 

where the numerical prefactor is lim n ^o(c/2/i2)- Our result equation 1171 then 
implies that V has the same dependence on its arguments as does the O(n) term 
in the CFT correlation function 

(T(2l)0l,2(23,23)</>l i2 (24,24)) . (35) 

This result generalises to other correlation functions, and implies that we may 
interpret the spin-2 component of the //-mass of loops which pass between two 
nearby points as being, in some sense, the derivative T of the CFT stress tensor 
with respect to n at n = 0. With this definition we then find the result 

(f ( Zl )f(z 3 )) = ' C ' i0)/ X ■ ( 36 ) 

as expected, where the left hand side is 

5 2 lim e- 4 J d8 12 e- 2 ^ J d6 M e~ 2 ^W ZlZ4 \ Z2Z3 (37) 

It would, of course, be important to establish this interpretation directly from 
the restriction property 



6 Twist operators in a simply connected domain 

The model may also be considered in a simply connected domain. All such 
domains may be conformally transformed to the upper half plane, with the real 
axis being the boundary. We may therefore derive the results for the upper 
half plane. The two point function of twist operators in the presence of this 
boundary satisfies the same set of partial differential equations as the four point 
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function in the bulk, with (23,23), (24, 24) assigned as the complex conjugates 
of the positions of the operators. This is a well known result in the theory of 
boundary conformal field theory [5] (BCFT) and follows from the condition of 
T = T on the boundary. The solution for the two point function is 

<to)^ 2 )>B CFT = ( 7 — ^; z ^r% — /8 " x 

M^i - z 2 ){z 1 - z 2 )(z 2 - z l ){z 1 - z 2 )J 



A( K ) [ Fr(l -j2-^;2-^) + B(n)(- V (l r?))^ 2 " 1 2 F t { J f - 1; \ ; r,) 

(38) 



where now there is only one cross-ratio 



'/ 



-B(/t) can be determined by looking at the boundary conditions as the two points 
approach the boundary, which is r/ — > — 00. Equation (|38[) may be analytically 
continued to large 77: 

( {z3 - Z ^- Zi) Y~^A(k)x (39) 

rzi^f wi K «.«.i, r r(2-f)r(i-f) r(|)r(i-§) 
V 211 4'4'2'Vlr(2-f)r(i-f) + r(f)r(i-f) 
r _ 3k/ 8 -i F /o _ 3« £-«. ..1 / r(2-f)r(§-i) r(f)r(f-i) 
( ?j 2 M 4' 4 > 2 2'ryH r(i - f )r(|) + ^r(f-i)r(f) 

(40) 

The two natural choices of B(k) are those which pick out one or the other 
conformal block as the operators approach the boundary. In the limit n — > 0, 
there are no loops in the loop gas picture so the correlation function should tend 
to unity for all 77. This is only possible if the first term in the square brackets 
of equation l|40|) is not present. A (non-unique) choice of B(k) which satisfies 
this requirement is that which picks out the second conformal block only: 

B{ r(2-f)r(i-f)r(f)r(i-f) 
r(f)r(i-f)r(2-^)r(i-f)' 

It is important to note also that this boundary condition is not compatible with 
the vanishing of the four point function in the limit 77 — > —00. The form of A(n) 
is dependent on the choice of normalisation of the twist operators. This can be 
seen from the r\ — > limit of equation (|38|l : 

km = AWf ^^ )^ 8 • 

A{k) must be of the form A = l + (rn + 0(n 2 ) due to the constraint that the four 
point function tends to unity as n — > 0. For small n, the correlation function 
may be expanded as 

(<Mzi)</>2(z2))bcft = 1 + n [t- ln <» - ^-^3^2(1,1,^; 2, |;»7) 

L37T 37T 3 3 

2Y(2\2 9 4 i 

+ (-»/(! - 2^1(3 , i; 3-5*7) + + 0(« 2 ) , 
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Gz!Z 2 () G(Zl)Z2 < ^Z 1 (Z 2 ) G( ZlZ2 ) 



{z\,Z\) (z 2 ,z 2 ) ((zi,zi) ) (z 2 ,z 2 ) (zi,z{) {(z 2 ,z 2 )) ((zi,zi) (z 2 ,z 2 ) 



(zi,zi) (z2,z 2 ) (zi,zi) (z 2 ,z 2 ) (zi,Zx) (z 2 ,z 2 ) (zi,z{) (z 2 ,z 2 ) 

Figure 7: The four types of configuration of a single loop around two points in 
the presence of a boundary. The complex conjugate points are included only 
for ease of calculation of quantities such as (— l)^ 1 ' 1 * (see text). 



where 

(zi - zf)(z2 - z 2 )a 



(zt - z 2 ){zl - z*)(z 2 - zl){z x - 4) ' 

In analogy with the interpretation of the four point function in the bulk in the 
loop gas picture described in section |4~T1 the following is the interpretation of 
the two point function in the presence of a boundary 

((-1)^.1* ) Ioopgag = (0(z 1 ,zT)0(z 2 ,2 2 -)) BCFT . (41) 

The left hand side is an expectation value in the ensemble of lattice loops, as 
before. N\^* is the number of times loops cross a defect line from z\ and z\. 
The set of all graphs may be decomposed into the set with no loops and the set 
with one loop; the other graphs are of order n 2 . The possible configurations of 
a single loop with a boundary are shown in figure [7] The two and one point 
functions are then found from the loop gas to be 

M = (0(z!)0(z 2 )) = ((-l)^' 1 * (-1)^-- ) = 1 - 2n(W 1{2) + W {1)2 ) + 0(n 2 ) 
Ci = {cj>{z x )) = ((-l)^ 1 *) = 1 - 2n(W {1)2 + W {12) ) + 0(n 2 ) 
C 2 = <</>(z 2 )) = ((-l) N ^-*) = 1 - 2n(W 1{2) + W (12) ) + 0(n 2 ) . 

Just as for the case of the loop gas in the bulk, the /z-masses {Wi} are the 
number of weighted loops belonging to the configuration i (sec figure . They 
are defined as 

w 1{2) = 5>i. 

Gl(2) 

Of the four possible configurations, the only one finite in the limit a — > is 
W {12) : 

W {m = ^^ + 0(n). 

In terms of r], 

W (12) = ~ lnfa(l ~ V)) - ^3^(1, 1, \; 2, §; rf) 

)2 -(- v (l- v ))i 2Fl (l,l4^) + 0(n). (42) 



67rr(4/3) v " "' * 1V 3' '3' 

Note, as for the bulk case, that the term involving a in M is cancelled by the 
subtraction of the product of one point functions. 
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7 Interpretation as a stochastic process 



We have seen in section 14.11 that the /i-mass of loops which wind around two 
of four points is invariant under conformal transformations. In this section, we 
show that the differential equations they satisfy can be interpreted in terms of 
an SLE process. Recall that these functions were identified as semi-connected 
four point functions, for example equation (|26|) 



: 4 \z 2 z 3 = lim(8n) 1 (<£(zi, zi)(f>{z 2 , z 2 )cj>{z 3 , z 3 )^>(z4,, z 4 )) 

n^O L 



C12C34 — C13C24 + C14C: 



■J:i 



where C^y are the correlation functions of a pair of twist operators at the points 
(zi,~Zi) and (zj,~Zj). The two point function is fixed by scale invariance and 
the four point function was determined as a solution to the partial differential 
equations 

r3,2 
2 21 



SJ and JTSJ. Using these, we find that W ZlZi \ Z2Z3 satisfies the PDE 



E 

Mi 
1 



1 



Z1Z4 \z2Z3 



1 

247 



Z2 - Zl \ Z2 - z 4 Z 3 - Z 2 



(z 4 

1 



Z 3 ~ Z 4 \Z 3 - z 4 



Z2 - Z3 



z 4 - z 1 \z 4 - z 3 



z 4 - z 2 



0(n) 



(43) 



together with a similar equation in which all the z s are replaced by z. Note that 
in the BPZ equations z and z can be taken as independent complex numbers, 
and it is only in applying these equations to physical quantities that one needs 
to impose reality conditions. These equations have almost the second order 
linear form which would result from applying the Ito formula to a martingale 
of an SLE process started at z\. One difference is that they arc complex. Wc 
can obtain a real equation by, for example, taking the real part of the sum of 
the holomorphic and antiholomorphic equations (the more general case will be 
discussed below), but note that when wc do this we have the sum 

d 2 Zl + dl = (d Zl + d- Zl f 2d Zl d- Zl = d 2 Xl ±A Z1 , 

where x\ = 3ie{z\) and A Zl is the Laplacian operator 9^ + dy ± . We also note 
that we can remove the inhomogenous part in H43(l by defining 

W = W / -(247r)- 1 25Re[-21n(z4-^i)+ln(23-Z4)+ln(z2-z 4 )-ln(z2-23)] (44) 
Then we can rewrite the real part of (|43|) as 



3^ i+ 25fte]T 



2d z 



Zl 



W = -A Z1 W. 



(45) 



Consider now a sequence of conformal maps gt{z) which satisfy the stochastic 
chordal Loewner equation 



dt 



9t{z) - z lt 



(46) 
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(Z2: Z 2 ) 



(23,23) 




1 Re 

Figure 8: A chordal SLE 8 / 3 growing in the half-plane ^sm(z) > 3m(zi) with its 
reflection in 9m(z) = ^sm(zi) 



where z\t = Z\ + \fHBt. Then if M(zi t , g t (zj)) were a martingale, its expec- 
tation value would satisfy (|45|l with the right-hand side set equal to zero. The 
coefficient of d\ would in general be k/2, so we should take k = 6. The right- 
hand side of i|45|) . together with the subtractions in (|44fl . therefore express the 
degree to which the weighted number of loops separating (211,24) from (2:2,2:3) 
fails to be a martingale. 

Eq. 021 describes a chordal SLE growing into the half-plane ^sm(z) > 3m(zi), 
together with its exact reflection in the line 3m(z) = ^sm(zi) as in figure[S] 11 . 
For k > 4 the hull, or outer perimeter, of this curve encloses a growing region 
of the plane. The entire interior of this hull at time t is mapped into the 
point i/JiBf. Because of this, we would not expect the weighted number of 
loops to be a martingale. Given some initial population of loops, each time 
that part of a loop is swallowed by the growing hull, it disappears from the 
population. The right-hand side of (|45|l must therefore express the rate at 
which this happens. It is reasonable that this should be proportional to A Zl W: 
if W were constant in some region, this would mean there were no loops passing 
through it. The leading non-zero rotationally invariant contribution should 
therefore be proportional to the Laplacian. 

Apart from this, the left-hand side expresses the fact that the /x-mass of loops 
which separate 212:4 and 2:223 and which have not yet been (partly or wholly) 
swallowed by the hull is the same as that in the conformally equivalent case 
of loops which separate (zu, 94(24)) from (g t (z2), (74(2:3)). This breaks down 
when the hull swallows any of the other three points. We conjecture that the 
subtractions in (|44|) take this into account. 

Of course, this is only suggestive and a number of important issues would have 
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to be resolved before one could actually derive our results from SLE. In partic- 
ular, one should explain why it is necessary to consider SLEg rather than some 
other value of k. This is presumably related to the fact that the hull of SLE 6 
corresponds locally to SLE 8 / 3 , and that both correspond to CFTs with central 
charge c = 0. The choice of a chordal SLE reflected in 3m(z) = 3m(zi) is 
clearly arbitrary; other choices correspond to taking different linear combina- 
tions of the holomorphic and anti-holomorphic equations. Perhaps using radial 
or whole-plane SLE would make the formula look more symmetrical. 

8 Conclusion 

In this paper we have studied scaling properties of loops in the loop gas picture 
of the O(n) model. In the picture, the partition function is a sum over all 
graphs of non-intersecting closed loops, weighted by nx l c where a; is a function 
of the reduced coupling and I is the length of the loop. We introduced twist 
operators, whose correlation functions count the loops separating the locations 
of the operators with weight n' different to the usual weight n, or equivalently 
count the minimum number of crossings of defect lines running between these 
locations. For the particular choice n' = —n, the twist operators have level two 
null states and their correlation functions satisfy BPZ type partial differential 
equations on the Riemann sphere. Thus, conformal field theory may be used to 
determine the analytic form of the two and four point functions. In the loop gas 
picture, the choice n' = ~n means that loops are counted are weighted by an 
additional factor of —1 to the power of the number of defect lines crossed, and 
the choice of path for the defect lines between the locations of the operators is 
unimportant. 

The limit n — + describes the theory of self-avoiding loops. In this limit, 
the partition function and correlation functions are dependent to first order 
in n only on the configurations of a single loop. Hence, by equating particular 
semi-connected parts of the four point function calculated using conformal field 
theory to the result from the Coulomb gas picture, we have deduced the expected 
number of weighted loops winding around two of the four points and shown that 
this number is invariant under conformal transformations. This is presumably 
the mass of such a subset under the measure on loops introduced by Werner. 
Other configurations of loops receive contributions from vanishingly small loops 
around a single point and are not finite in the limit of vanishing lattice spacing. 

The central charge of the O(n) model for small n may be found from the de- 
pendence of the /i-mass of loops around two of the four points as a function of 
the cross ratio of the positions of the four operators. The result c ~ 5n/3ir from 
the analytic results agrees with other methods of calculating the central charge 
of the O(n) model, and lends support to the interpretation of the stress tensor 
for curves satisfying conformal restriction given in Ref.JJ]. 

A similar calculation was also carried out for the model defined in a simply 
connected domain. A general simply connected domain may be mapped via 
a conformal transformation to the upper half plane with the real axis as the 
boundary. The two point function of operators in the upper half plane satisfies 
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the same partial differential equations as the four point function on the Riemann 
sphere, with two additional operators positioned at the complex conjugates of 
the original operators. By again equating the results from conformal field theory 
with those from the Coulomb gas picture, we have deduced the /x-mass of loops 
around the two points in the upper half plane. 

Finally, we have shown that the differential equations satisfied by the above 
quantities should have a stochastic interpretation in terms of a chordal SLEg 
process starting from one of the points Zy, along with its reflection in a fixed 
line. 
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